Doping induced spinless collective excitations of charge 2e in doped 

Mott insulators 
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Based on the Hubbard model description of the doped Mott insulator, we prove that there must ex- 
ist a spinless collective excitation of charge 2e and of energy U — 2(x peaked at momentum (n/a,n /a) 
in the non-superconducting state due to doping. Such a collective excitation arises from the electron 
pair correlation. Its existence is related to the breakdown of an approximate SU(2) particle-hole 
symmetry, from which an effective field theory of the new excitations is constructed. 



It has become clear that the normal state of the high 
T c cuprate superconductors shows many highly unusual 
properties. These properties are related to the fact that 
the high-r c cuprate superconductors are doped antifer- 
romagnetic (AF) Mott insulators. It is then not surpris- 
ing that the unusual behaviors are even more striking in 
the underdoped region, when the concentration of doped 
holes is small. A number of theoretical scenarios have 
emerged to understand the dynamics of the doped Mott 
insulators. Anderson jy introduced the notion of "spin- 
charge separation" in these materials and envisioned the 
resonating valance bond state — a liquid of spin sin- 
glet pairs. This notion of spin-charge separation natu- 
rally accounts for all the qualitative features of the spin 
gap (2). However, the spin-charge separation was proven 
only in one dimension whereas the cuprates are regarded 
as quasi-two dimensional. On the other hand, several 
scenarios of preformed pairs above T c || [| have been 
proposed as well to account for the underdoped (spin 
gap) regime of the high T c phase diagram. In this paper, 
we shall prove, based on the Hubbard model description 
of the high T c cuprates, that there must appear a novel 
kind of excitations that are spin singlet and carry charge 
2e in the non-superconducting state due to doping. Our 
results support the ideas of spin-charge separation and 
preformed pairs. An effective field theory of these charge 
excitations will be derived. 

Consider the extended Hubbard model in a 2D square 
lattice with spacing a, which we believe contains the es- 
sential physics of the doped Mott insulators. The Hamil- 
tonian is 

H = -t (cUx)ca(x 1 ) + h.c.) + J S(x) ■ S{x r ) 

(asas')fr (seas') 

+U^n ] {x)n i {x) - /i^n(x) , (1) 

X X 

where n a (x) = c\(x)c a (x) with spin a (=T, j), n(x) = 
ni(x) +ni(x), and /1 is the chemical potential that fixes 
the density of the system. 

The Hamiltonian (Q) is invariant under the transforma- 
tion, c„ — » Uaa'Ca'. This is the usual SU(2) spin symme- 
try. Furthermore, near half-filling the model has another 



(approximate) SU(2) particle-hole symmetry, sometimes 
called pseudospin group || ^, To see the symmetry, 
let us introduce two pseudospin doublets M 



ipi{x) = 



c-rO) 

(-)-c{(x) 



ip 2 (x) = 



ci(x) 
-(-) x c\(x)J 1 



where (— ) x = e lx Q with Q = (n/a,Tr/a). In terms of 
the doublets ip a , the Hamiltonian (|j) becomes 



(scsc') (assc') 



([7-2^)^03 



(2) 



where 



and 



<i> = -V'L'nmVw ■ 

The f is a vector of the usual Pauli matrices. (Note: 
Indices a, j3, • • ■ label two pseudospin doublets; l,m,- ■ ■ 
label the rows of each.) Explicitly, fa = \{n — 1), and 

A + , iA_ — ?A + ) where A_(ai) = 



6 2 ) = <^(A_ 



c±(x)c^(x) and A + = A I are the pairing fields. The 
theory is thus obviously invariant at half filling (/1 = 
U/2) under a global SU(2) pseudospin transformation: 
ip a l — » U/ m "0 Qm . Clearly, <j> transforms as a pseudospin 
vector. The operators 

J ± = Y.(fa± l fa) = J2(-) x A ±1 J 3 =X>, ( 3 ) 



generate the (approximate) SU(2) pseudospin symmetry, 
in the sense that 



[H,J±]=±{U-2/i)J ± , [H,J 3 }=0. 



(4) 



The symmetry becomes exact at half filling with fi = 
U/2 M. Note that J's commute with the ordinary spin 
operators S H. 
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Consider the non-superconducting, underdoped regime 
of the phase diagram. Hole doping breaks the SU(2) 
pseudospin symmetry down to the ordinary charge U(l) c 
subgroup in two ways. At first, it breaks the symme- 
try explicitly at the Hamiltonian level. This is obvious 
from Eq. (||). Secondly, doping forces the ground state to 
'line' up in the 3-direction (or z) in the pseudospin space 
while leaving the usual spin symmetry intact. This point 
is manifested by noticing that doping induces a ground 
state such that (<fo)o = — 7^ where Sh = 1 — («)o 
is the concentration of doped holes. The second way 
of breaking the symmetry is reminiscent of the notion 
of spontaneous symmetry breaking. One would imme- 
diately conclude from the Goldstone theorem jit)) that 
the Goldstone modes must appear at low energy. How- 
ever, the present case substantially differs from the case 
of spontaneously broken symmetries in that upon with- 
drawal of doping the ground state recovers the full SU(2) 
pseudospin invariance fll|| . The usual Goldstone theorem 
does not apply in the present case. We need to carefully 
examine the breakdown of the pseudospin symmetry. 

Collective excitations, if any, can be identified by 
examing the analytical properties of correlation func- 
tions. Consider the Fourier transform of the following 
(retarded) electron-pair correlation function 

D R + {x -x',t- 1') = -%6{t - 0([A- (sb, t), A + (x', t')}) . 

(5) 

At momentum k = Q, 

D\(Q,t- t') = -i6(t - t'){[J-{t), (-T'A + (x>, 0]>o , 

(6) 

where we have used the identity J_ = ^ x e t( 3' x A_(x). 
From (||), the time dependence of J_ can be determined 
via the equation of motion: 

J_(t) = e-W-^-t'U-it') . 

Also, we have the equal-time commutator 
[J-,(~) x ' A+(x')] = -2(f> 3 (x'). Therefore, we ob- 
tain exactly 

D* + (Q,u) = 2 fo>° , . (7) 

oj — (U — 2/i) + in 

The spectral function, defined in the Lchmann represen- 
tation as 

A(k,w)= (—) J2\(0\A-(0)\N)\ 2 6 2 (k-k N )6(u;-u N ), 

V a ' N 

_(8) 

can be determined exactly at momentum Q from (f7|): 

A(k = Q,uj) = -2ImL> R + (Q,w) 

= -4<^3)o7r5(o; - (17 - 2/x)) . (9) 

Thus as long as the symmetry is broken by doping with 
(^3)0 = — % 7^ 0, A{Q, to) cannot vanish, but rather con- 
sists entirely of a term proportional to 5(u> — (U — 2/i)). 



Such a term can obviously only arise in a theory that has 
particles of spectrum u> = U — 2/i at k = Q. Furthermore, 
a delta function S(u> — (U — 2/i)) can only arise from sin- 
gle particle states; multi-particle states would contribute 
a continuum. The operator A_ is bosonic, so (0|A_|iV) 
vanishes for any fermion state N. The state AL|0} is 
rotationally invariant in spin space, so (0|A_|A) must 
vanish for any state TV of non-zero spin. (We assumed a 
ground state of no long-range spin order; an AF order will 
change our conclusion.) Also (0| A_ |7V) vanishes for any 
state N that has different (unbroken) internal quantum 
numbers from A_. A similar analysis beginning with the 
correlation function D^_ in place of shall lead to 
another collective mode of energy — (U — 2/i) at k = Q. 
We then conclude that there must exist two (spinless) 
bosonic collective excitations of energy ±({7 — 2/t) at mo- 
mentum Q and of the same (unbroken) internal quantum 
numbers as A_ and A+, respectively. 

Let us denote two collective modes by W~ and W + , 
respectively, and consider how they should transform un- 
der the unbroken U(l) c symmetry. The quantum charge 
operator, if represented in terms of electron operators, is 
Q = — e Y^ x 71 = — 2e( J3 + M-) where M is the total num- 
ber of sites of the lattice and the electron charge is — e 
in our convention. Since W + (W~) must have the same 
U(l) c quantum number as A + (A_), the commutation 
relation of Q and C± has to be identical with that of Q 
and A±. Therefore, 

[Q,W ± ]=T2eW ± . (10) 

Physically, this means that W~ and W + create bosonic 
excitations of charges 2e and — 2e, respectively. Obvi- 
ously W + and W~ are conjugate, so there exists in fact 
only one kind of such modes. (Hereafter, we shall simply 
use W to denote the collective excitation.) We summa- 
rize these results into 

Theorem. If the ground state of the Hubbard model 
is such that (^3)0 = — % 7^ but (A±)o = 0, then there 
must exist a (spinless) collective excitation of charge — 2e 
and energy \U — 2/i|. 

We note that the above collective excitation may be 
viewed as a 'pseudo-Goldstone-like' mode from the break- 
down of the pseudospin symmetry. But we emphasize 
that the present case conceptually differs from that of 
the spontaneously broken symmetry, as discussed above 
(see also the note 0]). Since /t = {7/2 when half-filled, 
we expect that the energy of the collective excitation, 
\U — 2/i I, is low for small Sh- As we shall show below 
from the view-point of effective field theory, the energy 
of this collective mode disperses quadratically away from 
k = Q. 

We would like to comment that the (spinless) collec- 
tive excitation discussed in the present paper is related 
to the triplet of collective excitations of Zhang and 
the 7T excitations (though of spin quantum number 1) in 
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the unified S0(5) theory However, the triplet was 

proven to exist in the superconducting state only, and 
Zhang H argued that there is no direct way to couple 
to J + pair excitations experimentally unless the ground 
state is superconducting. By contast, our collective exci- 
tation can physically show up in the non- superconducting 
state with 5h acting as the role of order parameter. For 
instance, Eq. (0) suggests that the collective excitation 
here should be experimentally observable in the eletron 
pair response function at finite doping while the pairs 
need not to be condensed. This is the main result of the 
present paper. 

It is quite tempting to find an expression of the state 
of the collective excitation W in terms of electron states. 
The spectral function (||) peaking at momentum Q sug- 
gests that (A_(0)|W(fe)) / 0. So a naive way of look- 
ing for a state \W(k)) is perhaps out of the states 
A^_(x = 0)|0), projecting out all states of momenta dif- 
ferent than k. We can then write the state \W(k)) in a 
general form 



\W(k)) 



d 2 qf(q) C \(k-q)c\(q)\0) 



(11) 



BZ 



f(q) is some function of q undetermined. The terms 
indicated by '• • ■' represent all other possible states of 
appropriate quantum numbers. So a simple W boson 
may be viewed as a pair of electrons of opposite spins 
moving with a center-of-mass momentum k, and is not 
condensed in the normal state (i.e., (W)q = 0). This is 
consistent with the picture of preformed pairs in the non- 
superconducting state (|, ||]. The formation of bound 
states of holes has been seen numerically in a hole-doped 
antiferromagnet G3]. 

But one cannot simply identify the operator W with 
A_. By a direct caculation, 



[A_(x),Al{x') 



-2(j) 3 (x)S x ^ x i 
[S h - h'(x)]6 a 



(12) 



where h' = n — (n). Therefore, A_ cannot describe 
a true Boson since it fails yielding the standard Bosc- 
Einstein statistics, as required of W as a collective exci- 
tation. Rather, Eq. (12) implies that W should describe 
the collective motion of electron pairs with the quantum 
fluctuation h' integrated out. This point shall become 
a little more manifest in the language of effective field 
theory. 

An effective field theory of the W collective excitation 
can be derived from the breakdown of the pseudospin 
symmetry: SU(2)— >U(1) C with (03 )o ^ 0. According 
to the effective field theory approach of broken sym- 
metries pi ]) the (Goldstone-like) charge bosons may be 
identified (apart from normalization) with variables pa- 
rameterizing the coset space SU(2)/U(1) C . We can write 
ip a as an SU(2) transformation acting on ip a : 



where ip a is chosen such that 



1 7 + 



4>=-^ a rip a = (0,0, a). 

Accordingly, the is expressed as <fr a (x,t) = 
R a s(x, t)a(x, t) where R a f,(x, t) is a 3 x 3 orthogonal ma- 
trix determined by U _1 T a ll = R a bT~b- Clearly, a(x,t) — 
\/}2a 4> a (x, t) 2 . In place of the field variables ip a , our 
variables now are tj) a and whatever other variables are 
needed to parameterize the transformation \}{x,t). We 
choose 



U(x,t) = [1 - l J2r r (-rUx,t)}/\/ 1 + C(x,ty (14) 

r=l 

with 



R% = 0l/(03+<7). 



(15) 
(16) 



The lattice factor (— ) x is inherited from the definition 
of 4> r (see text following Eq. (^)) and is factorized out to 
ensure not to double the period of the lattice. The field 
variables £ r can be combined into a complex scalar field 

W = & + iC,2 , 

which may be identified as the collective excitation whose 
properties have been described in the theorem. 

As an example, consider a simple case in which 
holes are uniformly distributed: (a(x))o = (03(a;))o = 
— 4^. Also, we would generally expect a non- vanishing 
expectation value of the valance bond field Xxy = 
t{cl(x)ccr(y))o [15[- Starting with the Lagrangian cor- 
responding to (E), we replace the i/j a fields at each space- 
time point with Eq. ( |l3| ) and keep terms only relevant to 
the collective excitation field: 



W 



J w 



-2.V - W * idJ 

v {x* xy [l-W*(x)W(y)]+c.c.} 
^ ^[l + \W(xW}[l + \W(yW} 



+a(U-2fi)J2 



\W\ 2 - 1 
\W\ 2 + 1 



(17) 



ip a (x,t) = U(x,t)ijj a (x,t) 



(13) 



This gives the effective theory of lower energy excitations 
about the ground state we assumed. Substituting all a 
and x fields with their expectation values and expanding 
the Lagrangian in powers of W and W* , we have 

L cS = S h J2W*[id t - (U-2/j)]W- 2Re Xxy \W(x)\ 2 

X ( X y) 

+ Ei^ r («+"-] + -- (18) 

{xy) 
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The terms indicated by '• ■ •' will contain high powers of 
the W field (A detailed justification of the above effec- 
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field values of \ have been obtained for different phases of 
the Heisenberg- Hubbard model For a real uniform 
Xxy = Xi the Lagrangian (|l8|) yields the lowest order 



energy spectrum E{k) = |^[2 + cos(fc x a) + cos(k y a)] + 
(U — 2/i). Other phases of \ would yield different forms 
of energy spectrum but all should have energy {/ — 2/i at 
{ir/a, 7r/o). 

In conclusion, we predicted a spinless collective exci- 
tation of charge — 2e in the non-superconducting state of 
the underdoped high T c cuprates by analyzing the elec- 
tron pair correlation. Identifying these new excitations 
thus offers an opportunity of testing the validity of the 
Hubbard model for the high T c cuprates. Our result is 
consistent with the proposals of spin-charge separation |j[] 
and preformed pairs O, [| , which are widely used to ex- 
plain many important features of the normal state. 
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